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1. Introduction 

The Monte-Carlo simulation of lattice QCD with dynamical fermions is rather expensive. 
Therefore ideas and concepts are sometimes tested in low dimensional toy-models which are com- 
paratively cheap and often provide interesting insights. Two of these toy models are the two dimen- 
sional Gross-Neveu and the chiral Gross-Neveu model. They are renormalizable, asymptotically 
free and have a rich particle spectrum. The masses of the particles are generated dynamically. The 
chiral model is invariant under axial U(l) transformations, the standard Gross-Neveu model only 
under a discrete Z2-subgroup. Both models were treated in the large-Af approximation [[I], ||, |3j] and 
in perturbation theory [Q, ||] in the past and as both models are integrable also some exact results 
are available [0, 0j. Different discretizations were used to simulate the discrete model [§, ph. To 



our knowledge the chiral model has so far only been simulated with staggered fermions Q10| ] in 



3d. Here we present first results of simulations with Wilson fermions with which the restoration of 



chiral symmetry requires more care Qll| ] 



2. The chiral Gross-Neveu model 

The chiral Gross-Neveu model describes N flavors of self-interacting fermions in two dimen- 
sions. The standard way of writing its euclidian action is 

where the summations in flavor space are left implicit. In two dimensions such quartic interaction 
terms are renormalizable. Moreover the model is asymptotically free which makes it in a sense 
similar to QCD. 

An equivalent action is given by 

s = fd 2 x ^ + a + ir 5 n+4W+^^ + ^) ■ (2-2) 

Here real auxiliary fields a, IT and An have been introduced which are additional bosonic integra- 



tion variables in the path integral. If these are integrated out the action fl2.1| ) is recovered, hence 
both formulations are equivalent on the level of the generating functionals for fermionic correlation 
functions. The language with auxiliary fields is better accessible to numerical simulations because 
the fermionic fields enter bilinearly and can be integrated out. 

3. Lattice formulation with Wilson fermions 

On the lattice the action of the model with Wilson fermions takes the form 

S = a 2 Y,\Y[Dw + m + a + iY 5 n+4}Y+^2+^2+^ L ) • (3-1) 
x I *8s ^8p ^8v ) 



The Wilson operator is given by 



D w = ~ [Yn (V; + V M ) - aV; V M ] , (3.2) 
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with the forward and backward lattice differences 

1 



v*v(*) 



yf(x) 



\j/(x + ji) - \[r(x) 



-iaOji /L 



W{x-fL) 



(3.3) 
(3.4) 



We incorporate a phase factor into the definition which is equivalent to a particular choice of bound- 
ary conditions. For instance 6=0 corresponds to periodic and 6 = % to antiperiodic boundaries. 
The fermionic matrix D = Dw + mo + a + iysTl + 4 is neither hermitian nor anti-hermitian and 
unlike in QCD not even 75 -hermitian. For periodic and antiperiodic boundaries however it is real 
(in a Majorana representation of the y-matrices). 



We simulate the model with a Hybrid Monte-Carlo algorithm Q12Q. For this we introduce one 
complex pseudo fermion field (j) for each two flavors in order to represent the fermionic determi- 
nant as well as for each auxiliary field a conjugate momentum field. In the simulation of the two 
flavor theory most processor time is used for the integration of the equations of motion (leap frog 
integrator) and in particular in the solution of linear systems {DD^)x = (j) which is done with the 
conjugate gradient method. 



4. The Schrodinger functional 

In computer simulations finite volume is unavoidable. We can exploit it by using finite volume 
renormalization schemes, a particularly successful one in QCD being the Schrodinger functional 



scheme Q13) , [14J . We simulate the chiral Gross-Neveu model in a finite box with spatial extent L 
and temporal extent T. In the spatial dimension we apply (anti-)periodic boundary conditions and 
in the temporal one we have Dirichlet boundaries 

VA(0,xi)P_=p(xi) P + v(0,xi) =p(*i) 

y(T,xi)P+=p'{x 1 ) P-y(T,x l )=p'(x l ), (4.1) 

where P± = (1 ± 7o)/2. Such boundary conditions may in principle cause new divergences, which 
would require the introduction of boundary-counterterms into the action, but no such terms are 



necessary in our case []15[]. 

An observable is a polynomial in the functional derivatives with respect to sources f] , 77 and 
boundary fields 

dTj(x) 07](x) 

^ = -7^ ^ = j^n (4 - 3) 

dp(xi) op(xi) 

£'fr> = -*^rT ?W = xi£-^> (44) 

dp'(xi) op'(xi) 
which acts on the generating functional. 
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5. Chirally symmetric continuum limit 

In the formal continuum model a Ward identity associated with the axial U(l) symmetry can 
be derived 

-5— (A M (*)tf(y))=0, x + y, (5.1) 

where 6 is some arbitrary operator. 

In our lattice model chiral symmetry is broken explicitly by the Wilson term and we have to 
introduce a bare mass mo ^ and non-degenerate couplings g 2 P ^ g 2 s into the action (3.1) which 
increases the number of bare parameters to four: mo, g s , g 2 P and g\. Let us assume for a moment 
that we can define two different renormalized finite volume observables g \ and g2 and that we know 



two operators G\ and G%, which make ( |5.1| ) vanish only due to chiral symmetry and not due to some 
other symmetry of the action (e.g. parity or charge conjugation). We want to take a continuum limit 
of some other observables e.g. the step-scaling functions of g\ and g 2 at a fixed physical size of 
the lattice. We then need to tune the four bare parameters for a series of increasingly larger lattices 



in such a way that g\ and g2 are kept at some fixed values and at the same time (5.1) is imposed 
for both operators G\ and The continuum-extrapolated values of the observables would be 
universal predictions of the chirally invariant theory. 

At first sight the tuning of four parameters seems impracticable even in a two dimensional 



system. But there is a formal argumentation in the continuum [16] that might help us a lot. There 
it is established that the combination g\ — g 2 s /N does not renormalize and hence can be set to an 
arbitrary constant 1 

^ =C onst + ^. (5.2) 
Although we do not see how the formal argumentation can be reproduced on the lattice one may 



hope that (5.2) holds also there up to cutoff effects. For the tuning of gp and mo we can use large-Af 



results [Q| and perturbation theory [|TJ] as a first guess. 

6. Observables 

In our simulations we measure the following correlation functions 

.3 



/a(*o) =-7 £ (Y(x)YoY5¥(xK{yi)Y5^{zi)) (6.1) 
/p(*o) = -^ I (?(*)y 5 yA(*)C(yi)y 5 C(zi)) (6.2) 



x\ ,yi,zi 

h = -|z I <Wi)C(yi)> (") 

/4 = -^i I (ZVi)y&'?&)?<yi)i&'C{zi)) . (6.4) 

yi ,zi ,z[ 



We believe that this is the reason why the vector- vector interaction term in the action is left out in most large-Af 



calculations. 
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Figure 1: Left: The quotient m versus the bare coupling g 2 s on a lattice of size LxT=12xl3. The bare 
mass was set to mo = trie +0.2. The other couplings were gp = 10gv = gs- Right: The second derivative 
of the logarithm of the Schrodinger functional with respect to 9 versus the lattice size. The couplings and 
the bare mass are kept constant at mo — 0.023, gs = gp — gv = 0.1. 



The quotient 



m(xo) 



g>o/a(*o) 



(6.5) 



2/pOo) 

is independent of xq up to lattice artifacts and vanishes when the Ward identity ( |5.1[ ) is satisfied. On 
the left hand side of figure [j] this quotient is plotted and compared with perturbation theory. The 
bare mass was set to mo = mi +0.2, where ml is the 1-loop result for the critical mass. 



We can use and fa to define a finite volume running coupling 



g 2cx (fa)R-L 



(0) 
4 ' 



(fa)R=Z l -fa 



4,._ U {f: 



(0)n2 



(/2) 



2 U2 



(6.6) 



where and are the tree-level values. The last equality holds if we renormalize fa by 
requiring that it takes its tree-level value 



(fa)R = Zlfa=ff ] . 



(6.7) 



Other renormalized observables can be obtained from derivatives of the logarithm of the 
Schrodinger-Functional with respect to the angle B 



d 

dd 

d 2 



90 



logZ 
2 logZ 







(6.8) 



e=o 



+ 



-^Lv r (*)[i-'tt]v'(*+i)^ 

-a 4 r 

f(x)[l - n ] y (a- + 1 ) xjr(y) [1 + ri] v(y - i) } 



(6.9) 
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The first derivative vanishes due to the parity invariance of the action. The second one is finite and 
its magnitude strongly depends on the "physical" spatial volume of the system. The right hand side 
of figure |l| shows how this observable approaches zero when the lattice size grows at fixed bare 
parameters. 

At the moment we are still in the process of testing whether these (or other) observables will 
allow us to define a line of constant physics accurately enough. 

7. Conclusions 

To recover a chirally invariant Gross-Neveu theory from a lattice model with Wilson fermions 
requires the careful tuning of four bare parameters. We have introduced and tested several observ- 
ables that can be used to define a "line of constant physics". First simulation results are encour- 
aging, but whether some continuum extrapolated quantity can be calculated accurately enough to 
justify the effort in this model remains to be investigated. 
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